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The conversion of light into free electron-hole pairs constitutes the key process in the fields of 
photodetection and photovoltaics. The efficiency of this process depends on the competition of 
different relaxation pathways and can be greatly enhanced when photoexcited carriers do not lose 
energy as heat, but instead transfer their excess energy into the production of additional electron- 
hole pairs via carrier-carrier scattering processes. Here we use Optical pump - Terahertz probe 
measurements to show that in graphene the latter process is unprecedentedly efficient. Our results 
indicate that the production of secondary electron-hole pairs - carrier multiplication - dominates 
the ultrafast energy relaxation of photoexcited carriers, prevailing over optical phonon emission in 
a wide range of photon wavelengths. This singles out graphene as a promising material for highly 
efficient broadband extraction of light energy into electronic degrees of freedom, enabling a new 
class of high-efficiency optoelectronic and photovoltaic applications. 

In semiconductors, photoexcitation usually yields one electron-hole (e-h) pair per photon with the pair energy 
equal to the band gap. The excess energy, which is the photon energy minus the band gap, is typically lost as 
heat, i.e. transferred to lattice degrees of freedom. For many applications, however, it is highly desirable to identify 
materials in which a photoexcited carrier generates additional electron-hole pairs, leading to an enhanced quantum 
efficiency [l|. The unique properties of graphene, such as its gapless band structure, flat absorption spectrum Q and 
strong electron-electron interactions, make it a highly promising material for broadband photon-electron conversion 
Indeed, recent theoretical work has anticipated that multiple electron-hole pairs can be created by carrier-carrier 
scattering following a single absorbed photon in graphene (carrier multiplication) [H, H| . 

However, the experimental observation of carrier multiplication in graphene has remained elusive. The main 
challenge stems from the ultrafast character of the two main competing energy relaxation pathways: carrier-carrier 
scattering and optical phonon emission. In the former process the energy of photoexcited carriers remains in the 
electron system, being disbursed among secondary hot e-h pairs created in the Fermi sea, whereas in the phonon 
emission process the energy is lost to the lattice as heat. While recent experiments have shown that photoexcitation 
of graphene can generate hot carriers @, 0] , it remains unknown how efficient this process is with respect to optical 
phonon emission. Here, we employ an Optical pump - Terahertz (THz) probe measurement technique to elucidate the 
phenomena occurring in the energy relaxation cascade that follows after photoexcitation. In particular, we exploit 
the variation of the photon energy of the pump light to experimentally quantify the branching ratio between the two 
competing relaxation pathways. Changing the photon energy is crucial as it allows us to prepare the system with 
photoexcited carriers having a prescribed initial energy determined by the photon energy, and follow the ensuing 
energy relaxation dynamics. We show experimentally, in combination with theoretical modeling, that carrier-carrier 
scattering is the dominant relaxation process, providing clear evidence for highly efficient carrier multiplication in 
graphene. 

The employed technique consists of an ultrafast optical pump pulse that excites carriers and a THz probe pulse 
that passes through the sample after a variable delay time. The THz pulses afford an exquisite time-resolved probe 
of the high frequency photoconductive response of photoexcited carriers, as reviewed in Ref. @. This technique has 
been used before with a fixed pump wavelength to study charge dynamics in multilayer graphene (9j-[l3|. and with 
a variable pump wavelength to show that in some specific semiconductor bulk materials and nanostructures carrier 
multiplication can take place fl4Hl6| . The ability to vary the pump light wavelength provides a crucial advantage in 
probing carrier multiplication, since a higher pump photon energy will lead to an increased number of c-h pairs (per 
absorbed photon). Here we apply this technique to examine the energy relaxation cascade of photoexcited carriers 
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in monolayer graphene. We find that the photoexcited density of carriers with energy above the Fermi energy scales 
linearly with photon energy (for constant absorbed photon density). This scaling is found over a wide range that 
spans almost an order of magnitude in photon wavelength. Our experimentally observed linear scaling indicates 
that carrier multiplication is remarkably efficient, in excellent agreement with our results from a theoretical model 
that considers electron-electron scattering and electron-optical phonon scattering. We argue that for intrinsically 
doped graphene, carrier multiplication does not entail a redistribution of carriers between the valence and conduction 
bands as it does in conventional semiconductor systems 0, or undoped graphene [1, [f| ■ Instead, "hot carrier 

multiplication" occurs by elevating the carrier population above the Fermi level. These hot carriers can induce a 
strong optoelectronic response 0, @| ■ 

The graphene sample that we use for our study consists of a monolayer of CVD graphene transferred onto a 
quartz substrate. We use monolayer graphene, in contrast to previous optical pump - THz probe studies which used 
multilayer graphene |^4l3j]. From Raman spectroscopy we estimate a Fermi energy of fj, ~ 0.17 ± 0.05 eV, which 
corresponds to an intrinsic carrier concentration of ~ 2xl0 12 carriers/cm 2 . We further characterize the sample using 
THz transmission (without optical excitation) and find that the graphene monolayer has a spectrally flat absorption 
of ~ 5% in the 0.4-1.6 THz region. This absorption is due to intraband momentum scattering of the intrinsic carriers 
(Drudc conductivity) with an extracted average transport time of Tt r ~ 20 fs (see Appendix for a detailed sample 
characterization) . 

Our method for probing the energy relaxation dynamics of photoexcited carriers is illustrated in Fig. QJ,. We 
optically excite the graphene/quartz stack and examine the pump-induced change in THz transmission AT = T — To, 
were T and To are the transmission with and without photoexcitation, respectively. The THz pulses follow the pump 
pulses after a tunable time delay with a time resolution of ~100 fs. Because the graphene layer is thin, the change in 
transmission is directly related to the photoconductive response so that oc — Aer (see Appendix). Here Act is the 
pump-induced change in the THz conductivity, which we refer to as the THz photoconductivity. 

The essential features of the observed dynamics are exemplified by two typical differential transmission time traces 
shown in Fig. [TJd. The observed temporal evolution is non-monotonic: the differential transmission first rises in an 
approximately linear fashion, reaching a peak at a delay time T pea k ~ 200 fs, and finally decays on a longer timescale 
~ 1.4 ps. These observations can be explained as follows. During the initial rise, t < T poa k, carrier-carrier scattering 
between the photoexcited carriers and the carriers in the Fermi sea promotes carriers from below to above the Fermi 
level (see middle inset of Fig. QJ)). Thus the effective temperature of the carrier distribution increases, peaking at 
t ~ Tpeak when the photoexcited carriers have relaxed and a hot carrier distribution is established. Subsequent 
relaxation occurs due to electron-lattice cooling, which is a relatively slow process with a characteristic timescale 
(~ 1.4 ps) much greater than r pea k- 

The measured change in transmission ^r- is positive, meaning that the conductivity is reduced as a result 
of photoexcitation, Act < 0. The observation of negative photoconductivity is in agreement with recent THz 
pump-probe studies on monolayer graphene (l8l . [l9j . The reduction in conductivity is naturally related to secondary 
carrier excitation processes: Since the momentum scattering time T tr (e) increases with carrier energy e [20l - [22l ] . the 
creation of a hot carrier distribution after photoexcitation leads to a change in the real part of conductivity that has 
a negative sign (see Appendix). 

Key insight into the processes contributing to the energy relaxation cascade comes from examining how the 
differential transmission signal peak value scales with the photon energy hf ', which in our experiment is varied 
over a wide range, from the infrared (0.16 cV) to the ultraviolet (4.65 eV). Since the initial photoexcited carrier 
energy, = ^hf, determines where the cascade begins, we can track how the photoexcited carriers relax at each 
stage of the energy re-distribution process. As a first step, we analyze the peak value dependence on absorbed 
photon density iV p h ton (i-e. the number of absorbed photons per unit of area), which is displayed in Fig. [TJ:, for 
six different photon wavelengths. In the fluence regime employed here, the signal increases linearly with absorbed 
photon density for each photon energy. The linear dependence of THz photoconductivity on fluence at fixed photon 
energy indicates that in this regime each photoexcited carrier acts independently from the other photoexcited carriers. 

Proceeding with the analysis we observe that increasing the photon energy at a fixed absorbed photon density 
leads to a larger differential transmission signal at the peak. This is clear from the slopes in Fig. QJ: that increase 
with photon energy. The origin of the increased signal for increased photon energy (at fixed absorbed photon 
density) is shown schematically in Fig. [T}1. Here, increasing the photon energy leads to an increased number of 
electron-electron scattering events during the relaxation cascade (higher carrier multiplication factor) and thus a 
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hotter carrier distribution. In Fig. [2^, we show the effect of increasing the photon energy by plotting the peak 
differential transmission signal normalized by absorbed photon density. Notably, the normalized signal scales 
approximately linearly with the photon energy, whereas without carrier multiplication the normalized signal would 
be independent of photon energy. 

It is instructive to combine these observations in a unified picture which provides an intuitive "bird's eye view" of 
the energy relaxation cascade. We do this by plotting the interpolated experimental contours for the AT/Tq peak 
value as a function of photon energy and photon number (see Fig. HJl). Strikingly, the contours of constant AT /To 
bunch up at high photon energy and spread out at low photon energy. This confirms that the two ways to achieve a 
hotter distribution of carriers - either by increasing the absorbed photon density or by increasing the photon energy 
- are completely interchangeable. The hyperbolical shape also indicates that the differential transmission signal 
scales with the fluence (incident energy per area, iVphoton X hf). This constitutes a clear qualitative signature of 
carrier multiplication processes. Without carrier multiplication, the magnitude of the response would be determined 
only by the absorbed photon density, and not the excitation energy; the contour lines would have been essentially 
vertical, with no change in AT /To as photon energy is varied (see also Fig. [2f). From these "bird's eye view" plots 
we conclude that carrier multiplication plays an important role in the energy relaxation cascade. Below we develop 
this notion more quantitatively and estimate the (energy dependent) efficiency of carrier multiplication in graphene. 

The observed linear increase of the peak photoconductivity with hf is a clear signature of carrier multiplication. 
The efficiency of carrier multiplication depends on the branching ratio between the two processes: (i) electron-electron 
scattering vs. (ii) electron-optical phonon scattering that occur during the rise stage, < t < r pca k. To extract this 
branching ratio, we develop a simple model for energy relaxation in the photoexcitation cascade and compare it 
with the data. The relaxation can be described in a general form via de/dt = — J7ei- e i(e) — ^7oi-ph(e), where e is the 
photoexcited carrier energy and J c \- C \ and Jk-ph represent the energy relaxation rates for processes (i) and (ii). Our 
analysis relies on the fact that the characteristic time of the photoexcitation cascade starting at the photon energy 
e = hf/2 (at t = 0) and ending at the Fermi energy e w fi (at t = T pca k) is much longer than the carrier-carrier 
scattering time. Indeed, typical values T pea k ~ 200 fs measured in our experiment (see Fig. lb) are considerably 
longer than the reported values for carrier-carrier scattering times, which are well below 100 fs Oil, SI- Thus the 
electron subsystem during the cascade can be described using an effective electron temperature which is distinct 
from the lattice temperature. Thermal equilibration with the lattice in graphene is slow |2J-[26(, and in our sample 
takes ~ 1.4 ps (independent of photon energy and fluence in the regime considered here). This separation of time 
scales allows us to describe the electronic system using the electron temperature approximation. 

Carrier multiplication leads to the creation of secondary e-h pair excitations through electron-electron scattering. 
These secondary pairs give a negative contribution to the THz photoconductivity, since the momentum scattering 
time Tt r (e) increases with carrier energy e [20l - |22l [. Accounting for fast thermalization of the secondary carriers, the net 
change in the THz photoconductivity can be expressed as — 2aA r p hoton«7ei-ei, where iVphoton is the absorbed photon 
density, and the pre-factor a is estimated in the Appendix (the factor of 2 accounts for two carriers (electron and 
hole) produced per photon). Integrating over the cascade, we obtain the conductivity at the peak: 



In the case that electron-electron scattering processes dominate, J e y-e\ 3> Jci-ph, Eq. IA5I directly leads to linear 
scaling. In this case, since a is approximately energy independent (see Appendix), we have Act ~ — aNp^^hf , im- 
perfectly linear scaling. In a more realistic regime, when Sel—e\ ~ J7ci- P h, the THz photoconductivity in Eq. IA5I 
is sensitive to the branching ratio of the two processes, a(e) = J e \-e\/ Je\--ph- Using J e \- e \ from electron-electron 
scattering events described in Ref. fl7| and writing J7" i- P h = j(e — uJo)Q(t — — Ep) we obtain the branching 
ratios plotted in Fig. 2b. Here ujq is the optical phonon energy (0.2 eV), 9 is a step function and 7 (in (ps) 1 ) 
is the coupling constant between electrons and optical phonons, which we use as a fitting parameter to obtain the 
branching ratio. Since the branching ratio has a strong energy dependence (approximately oc 1/e at high energies), 
we use the scaling of Act with photon energy as a sensitive probe of the magnitude of the branching ratio. Indeed, 
Acr/iV p h from Eq. 1 begins to deviate from linear scaling (with hf) significantly when the branching ratio becomes 
smaller than unity, a(e) < 1. 

In Fig. 2a, we compare Eq. 1 for three values of 7 to our observed differential transmission signal normalized by 
absorbed photon density (see Appendix). The best-fit value 7 f» 4ps _1 corresponds to the solid curve in Fig. 2a. 
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The branching ratio values indicate that electron-electron scattering dominates the energy relaxation cascade. For 
stronger electron-phonon coupling, i.e. larger values of 7, we find that Act from Eq. 1 bends down at higher energies, 
significantly deviating from linearity, as illustrated by the dotted curve for 7 = 10ps _1 . We can exclude fits with 
7 > 10ps -1 as the data clearly lie in the range lps -1 < 7 ^ 10ps _1 . Inspecting the branching ratios that correspond 
to 7 = 1, 4, and 10ps _1 (see Fig. 2b) we conclude that the energy relaxation cascade that produces our observed 
AT /To have branching ratio a(e) ^ 1. Hence, electron-electron scattering dominates the energy relaxation cascade 
of photoexcited carriers, prevailing over electron-optical phonon scattering. 

We find good agreement between our experimental observations and values predicted by the model. First of all, 
the branching ratio we find in Fig. [2Jd is in good agreement with the theoretically predicted branching ratio. Using 
the known value of the electron-optical phonon deformation potential [24l |. we compute 7 = 1.36 ps _1 (see Appendix), 
which lies in the range of 7's obtained by fitting the data. In the Appendix we also show that the experimentally 
observed magnitude of the signal is in reasonable agreement with the theoretical result for a based on the effective 
temperature model. Furthermore, it is interesting to note that the rise time of the pump probe signal in Fig. \Tjp - 
corresponding to the time needed for the energy relaxation cascade - is longer for excitation with a higher photon 
energy. Although our experimental time resolution is only marginally smaller than the rise time, the variation with 
photon energy indicates that the observed behavior at t < T pca k is not limited by time resolution. The observed 
dependence is in qualitative agreement with the theoretical prediction of l ong er cascade times at higher photon 
energy, where more carrier-carrier scattering events occur during the cascade [171 ] . 

Electron-electron scattering dominated cascade and large branching ratios in graphene seen above have particularly 
strong ramifications for the quantum efficiency of graphene. Using the branching ratios plotted in Fig. 2b, we calculate 
the fraction of the photon energy that remains in the electronic system after the cascade (integrated e-c efficiency) as 

Via) = ~ -#TT* * = hf/2. (2) 

e * Jo a ( e ) + 1 

We plot r](ei) for 7 = 1 ps _1 — 10ps _1 in Fig. 2c showing that more than 50% of the photon energy remains in the 
electronic system even for photoexcitation energies as high as hf = 3 cV. This shows that carrier multiplication in 
graphene is highly efficient. We note that our model may overestimate the efficiency at low energies, as close to the 
Fermi surface energy relaxation may depend on pathways that were not included in the model: acoustic phonons, 
fiexural phonons and substrate surface phonons. However, since these processes only become important close to the 
Fermi surface we expect their impact on the total efficiency to be small. 

For most applications the relevant figure of merit is the integrated efficiency r\ as it describes the fraction of light 
energy that is passed to the electronic system. However it is also interesting to examine how many hot electrons are 
created from a single incident photon. Our model predicts a carrier multiplication factor that scales approximately 
linearly with photon energy, N ~ e;//i [Vi\ . Taking into account the extracted efficiency of 80% for excitation 



with 3 eV pump light, we find a carrier multiplication factor of 9 17J, which is well above the values reported for 
semiconductor systems 0, • 

Our study reveals efficient carrier multiplication in graphene, thereby resolving a long-standing question about the 
relative importance of electron-electron scattering versus emission of optical phonons in the energy relaxation cascade 
triggered by photoexcitation. Crucially, the transfer of energy from photoexcited carriers to electronic degrees of 
freedom in graphene is efficient over a wide range of frequencies (from the UV to the infrared), unlike conventional 
semiconductor systems where the frequency range is limited by the band gap. Furthermore, the efficiency of carrier 
multiplication is expected to be highly sensitive to the doping level [l?} , enabling effective manipulation of the energy 
cascade pathways. Thus graphene enables enhanced quantum efficiencies and tunable energy transfer over a wide 
spectral range. 
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III. FIGURE CAPTIONS 

Figure 1 

a) Experimental observation of carrier dynamics. Photocxcitation of electron-hole pairs in graphene monolayer by 
a pump pulse is followed by a Terahertz probe pulse that provides a measurement of the pump-induced change in 
THz electric field transmission AT = T — T (red curves). The time dynamics of the photoexcited carriers is studied 
by varying the time delay t between the pump pulse and the probe pulse. The change in transmission AT, which is 
directly proportional to the change in THz conductivity, is used to characterize the degree to which secondary e-h 
pairs are produced as the photoexcited electron cascades down to the Fermi level, b) Time resolved carrier dynamics 
for two different photon energies: the pump-induced change in THz transmission AT/Tq as a function of pump-probe 
delay time t for fixed absorbed photon density iVphoton ~ 10 12 cm~ 2 . The insets show a schematic representation 
of the carrier distribution N(e) before (t < 0) and long after (t 3> r pca k) pump excitation, and the hottest carrier 
distribution N*(e) directly after the energy relaxation cascade (t = r pca k). c) Scaling of the differential transmission 
signal ^r- peak values, obtained from the time traces as in B at t = r pca k for six photon energies as a function of 
absorbed photon density iV p hoton- The lines are linear fits to guide the eye, showing that the signal increases linearly 
with absorbed photon density. For a given absorbed photon density iV p h ton a higher photon energy hf leads to an 
increased signal, corresponding to a hotter carrier distribution, d) The effect of varying hf on the energy cascade 
illustrated for two photon energies, hf > hf. Photocxcitation creates a primary electron-hole pair, and triggers a 
cascade of carrier-carrier scattering steps, where energy is transferred to multiple secondary electron-hole pairs in the 
conduction band generating a hot carrier distribution. The number of secondary pairs increases with photon energy, 
leading to a hotter carrier distribution and a larger observed signal. 
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Figure 2 

a) Extraction of the branching ratio between e-e scattering and optical phonon emission from comparison of 
the experimental data and model: The pump-probe signal peak value normalized by absorbed photon density 
-^photon features approximately linear scaling with hf, indicating that electron-electron scattering dominates the 
energy relaxation cascade (see text). Deviation from linearity is accounted for in the model by including optical 
phonon emission with a coupling strength 7 (see Eq.l). The best-fit curve (7 = 4ps _1 , solid line) and the curves 
that marginally agree with the experimental data (7 = lps -1 and 7 = 10 ps" 1 , dashed line) are shown, b) The 
branching ratio for the carrier-carrier scattering and optical phonon emissions pathways as a function of initial carrier 
energy for the three coupling strengths shown in a. Larger-than-one values indicate that electron-electron scattering 
is the dominant energy-relaxation pathway, c) The integrated efficiency for the carrier-carrier scattering pathway as a 
function of carrier energy for the same three coupling strength values: best- fit (solid line), and upper and lower bounds. 
The estimated efficiencies are well above 50%, indicating that a large fraction of incident photon energy is transferred 
to the electronic system through efficient carrier multiplication. d,e,f) "Bird's eye view" of the differential signal vs. 
photon energy and absorbed photon number density for experimental data (d), model best- fit (e) and (simulated) 
phonon-cmission-dominated cascade (f). For efficient carrier multiplication, the differential transmission signal is 
expected to increase with both absorbed photon density iVphoton and with carrier energy hf /2. The contour lines of 
constant as a function of -/Vphoton and hf /2 obtained by interpolation of the data (d) indeed show this behavior. 
The good overall correspondence with the model for a best-fit 7 value (e) further corroborates the conclusion that 
e-e scattering dominates the cascade. In contrast, simulation of phonon-cmission dominated cascade (7 = 100 ps _1 , 
f) shows significant departure from the data. 
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IV. APPENDIX 
Al. OPTICAL PUMP - TERAHERTZ PROBE TECHNIQUE 

The optical pump - Terahertz probe setup is based on an ultrafast amplified laser system with a pulse duration 
of ^50 femtoseconds, a center wavelength of 800 nm and a repetition rate of 1 kHz. Part of the output is used to 
create the optical pump pulses and part is used to create the THz probe pulses. The optical pump pulses have a 
photon energy that we vary using non-linear optical conversion. To create 400 nm excitation pulses, we use frequency 
doubling in a /3-Barium Borate (BBO) nonlinear crystal (9 = 29.2°, 1 mm) and for the 267 nm pulses, we use the 
sum frequency of 400 nm and 800 nm pulses in a second BBO crystal (9 = 44.3°, 0.2 mm). For the excitation pulses 
in the infrared, we feed the 800 nm pulses into a home-built optical parametric amplifier (OPA), which we tune to 
cither produce signal pulses at 1400 nm or idler pulses at 2100 nm. Finally, for the 8 /im pulses, we use a AgGaS2 
nonlinear crystal (9 = 50°, 1.2 mm) for difference frequency mixing between 1455 nm idler pulses and 1778 nm idler 
pulses. The pump pulses pass through a diffuser to ensure spatially homogeneous excitation of the graphene sample 
and through a 500 Hz chopper so that we alternately probe the THz electric field transmission with and without 
photo-excitation by pump pulses, To + AT and T , respectively. We determine the pump fluence at each wavelength 
by measuring the transmission trough five holes of different well-known sizes with a calibrated power meter. The 
8 fan and 267 nm pump pulses had very low fluences and therefore did not pass through the diffuser. Excitation 
without diffuser could lead to an underestimation of the absorbed photon density in the case of hot spots in the beam 
profile and therefore these two photon energies are not taken into account for the determination of the efficiency of 
carrier multiplication. 

The THz probe pulses are created through optical rectification in ZnTe nonlinear crystal (110 orientation, 0.5 mm) 
and focused on the sample. After recollecting and collimating, the THz beam is focused on a second ZnTe crystal and 
overlapped in time with an 800 nm sampling beam, whose polarization is rotated by the presence of the electric field of 
the THz beam. We then detect this polarization rotation, which is directly proportional to the electric field strength 
of the THz beam. By varying the time delay between the THz beam and the 800 nm sampling beam, we obtain the 
time evolution of the THz pulse. For the pump-probe measurements described here, we position this time delay such 
that we are at the peak of the THz pulse. We then change the time delay of the pump pulses with respect to the peak 
of the THz pulse in order to follow the temporal evolution of the THz pump-probe signal, which is proportional to the 
THz photoconductivity in the thin film approximation: Act = . Here n s = 1.95 [27] is the refractive index 

of the quartz substrate and Zq = 377 f2 is the free space impedance. All measurements were done at room temperature. 



A2. SAMPLE CHARACTERIZATION 

The graphene samples consist of a large (~ 1 cm 2 ) single graphene layer that was grown by chemical vapor 
deposition (CVD) and subsequently transferred onto a 200 /im thick fused silica substrate (CFS-2520 from UQG 
Optics). To determine the intrinsic doping of the graphene layer, we performed Raman spectroscopy (see Fig. I Alb ). 
We find a G peak located at 1586 cm -1 corresponding to a Fermi level of Ep = 0.17 ± 0.05 eV and thus an intrinsic 
doping of riintr ~ 2 • 10 12 carriers/cm -2 [3(| Hl|. The FWHM of the 2D peak is 30±5 cm" 1 , confirming that our 
sample is monolayer graphene. We measured the THz transmission (without optical pump excitation) and find that 
the graphene layer leads to a spectrally flat absorption of ~ 5% in the THz region (see Fig. lAlb ). Using the thin film 
approximation we find that this absorption corresponds to an intrinsic sheet conductivity of CTi, ltr « Ts ^ b ~ T ° = 
4-10" 4 S/m (or ~ 10e 2 //i), with T su b the transmission through the substrate without graphene. The (Drude) sheet 
conductivity in graphene is determined by the carrier density and the weighted momentum scattering time (rt r ). We 
use the carrier density from the Raman measurements to extract a momentum scattering time of ~20 fs, similar 
to earlier conductivity studies of large area CVD graphene [1, H)]. This shows that the extracted intrinsic carrier 
density is realistic and that THz radiation is an excellent probe for the free carrier conductivity. We note that the 
substrate thickness of 200 /im results in a small additional peak in the pump-probe delay traces at a delay time 
of ~3 ps due to reflection of the THz wave. We have confirmed that this peak is absent when using a thicker substrate. 

In our measurements we measure the pump-induced change of the conductivity of the free carriers, i.e. the THz 
photoconductivity. In order to make quantitative comparisons between the THz photoconductivity induced by 
photo-excitation with different photon energies, we need to accurately know the number of absorbed photons at each 
pump wavelength. We obtain this number through the measurement of the fluence (the energy of the pump beam 
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per unit of area) and through the graphene absorption at each pump wavelength. We determined the graphene 
absorption using a standard double beam spectrometer (Perkin Elmer Lambda 950), where we placed an identical 
quartz substrate without graphene in the reference path. The absorption was determined from 200 nm up to 2.2 /iin 
(see Fig. lAlb ). above which the substrate starts absorbing strongly. The change in reflection due to the presence of 
graphene is ne eligibly small (<0.1%) 0. We find relatively flat absorption that is slightly lower than the universal 
value of 2.3% 3, down to a wavelength of ~500 nm and increased absorption due to excitonic effects that peak 
around 250 nm 12911 . 



A3. DIFFERENTIAL TRANSMISSION, THZ CONDUCTIVITY, AND THZ PHOTOCONDUCTIVITY 

The transmission, T(lu), of light at a frequency u through a single gra phene layer and substrate stack (in the thin 
film approximation) is related to its optical conductivity, a(w) via [9H32j 



nub 



cr(w)Zo 



n s + 1 



(Al) 



where T su b is the transmission through the substrate only (no graphene), Zq = -\/eo/Vo is the impedance of free 
space, and n s is the index of refraction of the substrate. Throughout this text, we only consider the real part of the 
optical conductivity which we will denote by a. In pump-probe experiments, It is natural to consider the change in 
transmission that is due to the pumping of the sample. The difference between pump on and pump off transmission, 
AT = T — Tq where To corresponds to the transmission with pump off, is also related to the change in optical 
conductivity of the graphene sample due to optical pumping. This gives 

AT = l + W^ol/K + l] _ , M _ Mw) _^ > (A2) 
T 1 + [<7i(w)Z ]/[n a + 1] 1 V + l' V ; 

where ctq(uj) is the optical conductivity with pump off, Act = <j\{uj) — <jq{u>) and in the last line we have used 
[<7(u;)Zo]/[n s + 1] < 1 and expanded the denominator. o\ is the conductivity when the pump is on. For uj in the 
ranges of 1 THz this is referred to as the Terahertz conductivity and probes intra-band scattering events as fiui <C Ep- 
The optical conductivity depends directly on the energy- dependent carrier distribution and gives a different value 
according to how far away from equilibrium the carrier distribution is. To see this, we use the kinetic equation to 
derive the conductivity [33j 

<9 t n k + eEV k n k = W k ',k [(1 - n k )n k ' - (1 - n k /)n k ] = ^ W k ', k (n k / - n k ) (A3) 

k' k' 

The perturbation of the non-equilibrium distribution due to an external electric field, ?i k = n k + <5n k , can be found 
by the relaxation time approximation where we linearize the right hand side in (5n k giving PVk^k (n k < — n k ) ss 
— (5n k /T tr (k) where Tt ^ = J^ k , Wk',k(l — cos6>). At finite frequency the carrier distribution response to an AC field 

E oc e~ luJt gives a (real part of) conductivity 

k 1 + lwT tr (ek)J 2 

where we have used the formula for current j = ^ k ei' k <5n k and picked out the w response. 

As described in the main text upon pumping by optical pulse, a narrow band of photo-excited carriers is created. 
These photo-excited carriers can relax energetically in two primary ways depicted in Fig. 1 of the main text : (i) 
carrier-carrier scattering with carriers in the Fermi sea or by emission of optical phonons. Both processes can lead to 
a change in the THz conductivity : termed THz photoconductivity. The former excites the carrier distribution about 
the Fermi surface leading to a hotter carrier distribution and a change in the THz conductivity. Because the transport 
time increases with carrier energy (2fj| - [22l |. the change in the real part of conductivity (measured in transmission) is 
negative. The latter process produces optical phonons which can act as scattering centers for carriers in the Fermi 
sea, increasing the resistivity of the sample. 

In the following we consider both processes. Accounting for fast thermalization of the secondary carriers, the net 

change in the THz photoconductivity can be expressed as — 2iV p hoton^J?ci-ci + &»7ei-ph^ , with the energy relaxation 
of the photoexcited carriers described in a general form via de/dt = — ^J e \- e \(e) — J r e i_ p h(e), where e is the photo- 
excited carrier energy and J c \- C \ and J^i-ph represent the energy relaxation rates for processes (i) and (ii). Here a is 
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the coefficient governing the change in THz conductivity from the capturing of energy by the electronic system (see 
below), and b is the coefficient governing the change in THz conductivity from the emission of optical phonons that 
act as scattering centers (see below). As a result, the THz photoconductivity from these two processes is 

Act = -2N photon / dt{aj el - el + bj cl _ ph ) = -2N photon / dE cl ph , (A5) 

To get a gauge of which processes contribute most to the photoconductivity, we note that for a photon energy 
of hf = 3.1 eV, and iVphoton = 10 cm our experiments found a transmission change of 0.75% (Fig. lc of main 
text) which corresponds to a photoconductivity of Act ps — 1.5e 2 /h. Using values of b, Ep = 0.17eV in our samples, 
t/ei-ei, and J7ei- P h calculated below, we find that the optical phonon scattering contribution to photoconductivity is 
Affph ~ —0.02 e 2 /h about 75 times smaller than what is observed. Hence, we conclude that the scattering off emitted 
optical phonons produce negligibly small contribution to the photoconductivity. As a result, we set b = in our 
analysis of our data (see Eq. 1, Fig. 2 of main text and discussion below). 



A4. ELECTRONIC TEMPERATURE MODEL FOR a 



We consider the photoconductivity due to carrier-photoexcited carrier scattering. This process excites the carrier 
distribution about the Fermi surface. Because carrier-carrier scattering is fast on the order of 10s of fs [13, [H| - 
shorter than the initial rise time ~ 200 fs measured in our pump-probe experiment (see Fig. lc of the main text) - the 
excited carrier distribution thermalizes quickly and can be described by an increased effective electronic temperature, 
T. Thermal equilibration with the lattice in graphene is slow [2^, [26|, [34j , and in our samples takes ~ 1.4 ps (see Fig. 
lb of main text). This separation of time scales allows us to treat the electronic system as out of equilibrium with 
the lattice. 

In the degenerate limit, E F 3> k B T, we can use the Sommerfeld expansion in Eq. IA4I giving a (real part of) THz 
conductivity 



ct(c) a ctHt=o + ^v{E F )k 2 B T 2 ^^ 

Ot z 



F(e) 



2 2 

e v 



■=E f 



1 



nr(e) 



■ w 



[rtr(e)] 



(A6) 



where Ep is the Fermi energy , i/(e) is the total density of states in graphene, and v is the Fermi velocity. We have 
kept the total carrier density, n, constant by accounting for changes in chemical potential as a function of temperature 
k 2 B T 2 1 Ep. Here we have neglected the temperature dependence of the transport time, Tt r (e), since we 

6] and can be 



/' 



E f 



6 B 



estimate that fc^AT < Ep. We note that these, in principle, can provide additional terms to Eq. 
included in a more sophisticated analysis. 

Therefore, the hot carrier contribution to the photoconductivity arising from an increase in temperature of the 
carrier-distribution, Act c i = o\ — ctq is 



Aov, 



h 2 T 2 

K B 1 l 



h 2 T 2 

K B 1 Q 



(A7) 



where 1, subscripts indicate pump on and off respectively. 

We note parenthetically that the differential transmission in our set-up is related to the the photoconductivity above 

but with a sig n that is opposite to 9 ^ e i e ' . For graphene, d ^Jjf 1 < for small probing frequencies u> ~ THz because 
Tt r oc e 20h22 



As a result, our observed positive differential transmission (or equivalcntly negative photoconductivity) 
is consistent with a hotter carrier distribution. On the other hand, if carrier density had increased (for example, by 
exciting electrons from the valence band), this would have resulted in a negative photoconductivity (which we did not 
observe) . 

We can relate Act c i to the change in heat in the system. In the degenerate limit, the electronic heat capacity is 



C = aT where a = ^v{Ep)k 2 3 . Hence, the photoconductivity is 



Act, 



k 2 AQ ir 2 



v{E F ) 



3 2 F{e) 



AQ 



a 



(T 2 - T 2 ). 



(A8) 



where AQ = Qi — Qo- This means that Act i is a direct measure of the amount of heat captured by the electronic 
system. Since the amount of heat absorbed by the carriers in the Fermi sea comes directly from the energy-loss of the 
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photo-excited carrier, the rate of heat entering the Fermi sea is 

dAQ 



dt 



i/el-el-^photon (A9) 



where — J e i— e \ is the energy-loss rate of the photo-excited carrier due to carrier-carrier scattering between the photo- 
excited carrier and the carriers in the Fermi sea and iVphoton is the absorbed photon density (ie. density of initial 
photo-excited carriers) . We have neglected thermal equilibration with the lattice as the time scales for cooling to the 
lattice ~ 2 ps are far longer than the time scales of energy-loss of the photo-excited carrier ~ 40 fs [T3, HH . 

As a result, the photoconductivity from the relaxation of the photo-excited carriers produces the first term in Eq. 
IA5lwith 

— ^ (A10, 

Here a is a positive quantity. As seen above, the photo-conductivity (and hence the differential transmission) is directly 
sensitive to the energy relaxation dynamics of the photo-excited carrier and Acr e i a unique probe of the amount of 
energy that gets captured by the electronic system. 

In the following, we estimate a from this model. Since conductivity in graphene for high doping depends linearly 
on density n ot/i 2 [2^ |. the einstein relation, a = e 2 v 2 v(i±)T{n) /2 tells us that r(e) = ae, where a is a proportionality 
constant. The factor of 2 in the einstein relation comes from dimensionality. This is consistent with carriers scattering 
off Coulomb disorder in graphene [20Tj22j . We measured a transport time at e = [i of 20 fs (see above section) and 
a doping corresponding to Ef = 170 meV (see Raman Spectroscopy above). Therefore, we infer that in our sample 
a « 117fseV _1 . 

Differentiating F, and using Ep = 0.17 eV and a typical THz frequency that was used u = 27rTHz (corresponding 
to / ~ 1 THz) we obtain 



a k 2.68 x 10~ 12 x ^-cm 2 eV _1 (All) 

This value for a gives a photoconductivity of Act = — 0.83e 2 /h for a photon energy of 3.1 eV and absorbed photon 
density of iVphoton = 10 11 cm , using the theoretical values of J c \- C \ and J7ci-ph obtained below. Experimentally, 
we measure a differential transmission signal of AT/Tq ~ 0.75% for these conditions, which corresponds to Act = 
— 1.5e 2 //i, in very reasonable agreement with the theoretical result. A comparison of a with b (estimated below) yields 
a/b ~ 5 for photoexcited carrier energy of e = 1 eV and ji = 0.17 eV. 



A5. OPTICAL PHONON SCATTERING CONTRIBUTION TO PHOTOCONDUCTIVITY, b 



We now consider the contribution to photoconductivity that comes from the emission of optical phonons. Optical 
phonons can act as scattering centers for carriers in the Fermi sea and alter the mobility of the graphene sample. 
Hence, a larger population of optical phonons increases the resistivity of the sample, Ap p h, and can lead to a negative 
photoconductivity, Acr p h = -OnAppi,. Here, <tq is the conductivity without pump excitation. In particular, the optical 
phonon limited resistivity is [35j 

A Pph = 2 ° p 2 1 ° ! (A12) 

where D op is an effective electron-optical phonon coupling constant }35l ]. p m is the mass density of graphene, u>o is 
the optical phonon energy of graphene and N(ojq) is the occupancy of optical phonons in graphene. 

To estimate the occupancy N(ujq we note that the optical phonons that are emitted from the energy-relaxation of 
the photo-excited carrier similarly only have momenta < e/(hv) where e is the energy of the photo-excited carrier. 
Hence, the (conditional) probability density, p, of emission of an optical phonon with wavevector q given that an 
optical phonon was emitted by the energy relaxation of a high energy photo-excited carrier with wavevector k is 

d9 

d P = —{i- cose) = 3^ ^;;;; , ™ dq, | q | = 2|k| C ose 




where 9 is the angle between q and — k, and the photo-excited carrier's energy is e = u7i.|k|. Here we have used the 
factor 1 — cos# to account for the coherence factor in the scattering matrix element. Since optical phonon scattering 
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with carriers in the Fermi sea is mainly limited kp far smaller than the Brillouin zone we are interested in the fraction 
of emitted optical phonons that lie within the circle of radius 2kp. The fraction of optical phonons emitted that lie 
within the circle of radius 2kp is 

As a result over the course of the relaxation of photo-excited carriers, the number density of optical phonons emitted 
within this circle of radius 2kp, is 

£ N(\ q \) = N photon [ ■ 7cl - ph/ "° xlfefc (A15) 

|q|<2fcF 

where iV(|q|) is the occupancy of optical phonons emitted. On average, each optical phonon in this circle of radius 2kp 
contributes the same amount to the transport scattering rate between optical phonons and carriers in the Fermi surface. 
Hence, the average occupancy within this circle, kp = [J2\q\<2k F ^(lll)] /(J2\q\<2k F )' ai l° ws us to estimate the 

change in resistivity coming from these optical phonons via Eq. IA12I Therefore, noting that X)| q |<2/c F = ^(2^ = 
kp /tt and using Eq. IA15I and Eq. IA121 we find that emission of optical phonons contributes a photoconductivity 
described by the second term of Eq. IA5l with 



27ra 2 h 2 D 2 op 
e 2 p m vu!^kpe 



(A16) 



We can estimate the value of b. Taking D op = 22.4cVl 1 [II], p m = 7.6 x 10~ n kg cm" 2 , uj q = 0.2cV, and 
/j = 0.17eV we obtain 

5.74 xlO- 13 e 2 2 ! 

b ~ ; f-r3f7-— -7 . p t n \ x -rem cV , (A17) 

(/x[eV]/0.2) x (e[eV]) h v ' 

where we have used the measured conductivity ctq ~ 10(e 2 //i) (see section [A"2j) . 



This value for b gives a photoconductivity of Atr = — 0.02e 2 //i for a photon energy of 3.1 eV and absorbed photon 
density of JV p hoton = 10 11 cm" 2 . Experimentally, we found Act = —1.5e 2 /h, which is significantly higher than the 
signal due to phonon emission only. We showed in Section IA4I that the signal due to hot carriers is much larger and 
closer to the experimentally observed photoconductivity signal. 



A6. ENERGY-RELAXATION OF PHOTO-EXCITED CARRIERS, -J 

The energy relaxation of a photo-excited carrier occurs via two processes: (i) scattering between the photo-excited 
carrier and the carriers in the Fermi sea (see Fig. lb main text) allows the energy to be transferred to the Fermi sea 
to create a hot carrier distribution, and (ii) the emission of optical phonons. Because of the large number of carriers 
in doped graphene, the former process can be very efficient [171 ]. 

Additionally, the energy relaxation via this channel occurs in relatively small steps of order fi the energy relaxation 
rate from impact excitation events, Jel-el [I3]i i s 



Jd-ci(e) = (piieV}) 2 ^/^ eV/ps (A18) 

where £(e, /z) is plotted in Fig. lA2l using a numerical integration detailed in Ref. [l7| (starred points in Fig. IA2[) . We 
subsequently used a 4th order polynomial to interpolate between the points (dashed line in Fig. IA2[) . For our doping, 
fx = 0.17 eV we find an efficient relaxation rate that for typical carrier energies are a couple of eV/ps indicating 
that the relaxation of photoexcited carriers with typical energy of e = 1 eV via J^ci-ci occurs with a few hundred 
femtoseconds. This is consistent with our experimental observation of a rise time of around T pea k ~ 200 fs (see main 
text Fig. 1). 

An alternative channel for energy relaxation of the photo-excited carrier occurs through the emission of optical 
phonons and gives an energy relaxation rate of — jT e \-ph- The transition rate of this process [36| can be described by 
Fermi's golden rule 
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K7 h = — E l M ( k '' k )l 2(5 ( Ae k',k + w q )^, k+q (iVK) + 1) (A19) 
q 

where Ae k ',k = £k' — £ki = = 200 meV is the optical phonon dispersion relation, and 7V(w q ) is a Bose function. 
Here k is the initial momentum of the photo-excited electron, k' is the momentum it gets scattered into, and q is the 
momentum of the optical phonon. The electron-phonon matrix element M(k',k) [3~H is 

|M(k',k)| 2 =. 9o 2 F k , k ,, g = ^ = (A20) 

where i^k.k' is the coherence factor for graphene, go is the electron-optical phonon coupling constant, a is the 
distance between nearest neighbor carbon atoms, and p m is the mass density of graphene. The energy-loss rate of the 
photo-excited carrier at energy e due to the emission of an optical phonon is 

^el-ph(e) = J2 <7k Ph (4 - e) [1 - /M] (A21) 

k' 

Integrating over q and k' we obtain 

Jei-ph(e) = ^Y^gl [l - f(e - wo)] (N(oj ) + l)u{e - oj ) = 7 (e - Lj Q )6(e - Ep - wo), (A22) 

where v[e) = e/ (2-kv 2 Ti 2 ) is the electron density of states in graphene and we have approximated (N{loq) + 1) « 1 and 
1 — f(e — luq) m 0(e — luq — Ep). 7 is a constant determined by the electron-optical phonon coupling. Hence, J7 e i_ p h(e) 
varies linearly with the photo-excited carrier energy e > w and vanishes for e < cj . Because the electron-phonon 
coupling with optical phonon is a constant, this result is to be expected from the increased phase space to scatter 
into at higher photo-excited carrier energy. 

Using p m = 7.6 x 10~ n kg cm" 2 and a = 1.42A in Eq. IA20I above yields 7 ps 1.36ps _1 . 



A7. COMPARISON OF DATA WITH MODEL 



These theoretical results for Je\- e i and Je\-ph lead to a branching ratio between electron-electron scattering 
and electron-optical phonon scattering that we wish to compare to the experimental data. Therefore we take 
the experimental results for the signal normalized by absorbed photon density (Fig. 2a of the main text), which 
is directly proportional to — -j^, and compare these with the model through Eq. 1 of the main text. We hx 

j7ei-ei to the theoretical result and use J7ci-ph = 70(e — Ep — ojq), where we use the electron-phonon coupling 
constant 7 as the parameter that determines the experimental branching ratio. The parameter 7 reflects how 
— -j^ 2 - depends on carrier energy (or hf): A low value of 7 corresponds to an energy cascade that is dominated by 

^ ' ph 

el-el scattering and results in linear scaling with carrier energy, whereas a high value corresponds to el-ph domi- 
nating, which leads to the signal bending down; large 7 leads to a non-linearity in Act setting in at a low carrier energy. 

We repeat the fit routine for a range of values for 7 and include one fit parameter that determines the amplitude 
of the signal, thus serving as a mere scaling factor. In Fig. 2a we use the amplitude obtained from the best 
fit for all three curves, such that all curves for different 7 pass through the same data point at low excitation 
energy (0.56 eV). We do not include the lowest (0.16 eV) and highest (4.65 eV) photon energies. In the former, 
the model does not cover this region and the absorbed fiuence is not known very precisely. In the latter, the 
diffuscr could not be used, which leads to possible over-estimation of the absorbed photon density (in case of 
hot spots in the beam profile). The fit results for 7 = 1, 4 and 10 ps _1 are shown in Fig. 2a of the main 
text. Here we show the residuals (squared difference between experimental data points and model) for a range 
of 7's, clearly showing a minimum around 4 ps _1 and strong deviation especially at large values of 7 (see Fig. 
S2 b). This clearly shows that our data is incompatible with a phonon-dominated energy cascade. The best 
fitting value is somewhat larger than the theoretically obtained one (1.36 ps" 1 ), but the theoretical value still lies 
well within reach of the experimental data. It is also conceivable that the data leads to an overestimation of 7 
(and thus an underestimation of the contribution of el-el scattering) due to a nonlinear (saturation) response of 
the signal at high photon fluences and high photon energy. This means that we extract a conservative branching ratio. 
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A8. APPENDIX: FIGURES AND CAPTIONS 
A | ! 1 
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FIG. Al: Sample characterization, a) Raman spectroscopy to determine the intrinsic carrier density through substrate doping, 
b) Terahertz transmission spectrum, showing spectrally flat absorption of ~5%, due to the conductivity of the intrinsic free 
carriers, c) Absorption spectrum to determine the absorption for the used pump wavelengths, indicated by black crosses. 
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FIG. A2: a) (Left ordinate) Starred points denote numerically evaluation of scaling function £(e//i) using Ref. [13]. Dashed 
lines are a 4th order polynomial fit of the numerical evaluation within this range : = -3.1081 x 10" 3 s 4 + 0.24 031a 3 - 

6.8417k 2 + 92.622a; - 59.383 where x = e/fx. (Right ordinate) Magnitude of J" e i_ e i for [i = 0.2 eV using Eq. |Xl8l b) The 
squared residuals as a function of 7, obtained by fitting the theoretical results for (e) to the experimentally obtained 

(e) oc -^-(e). The fit is done for a range of values for 7, which reflects the branching ratio, and with one free fitting 

scaling factor. The best fit is obtained for 7 =4 ps -1 , and the residuals increase strongly for large 



parameter that is just 
values of 7. 



